G3 Genes|Genomes|Genetics

Reporters are widely used in Drosophila genetics to visualize gene expression and cell lineages. However, uncharacterized limitations in specific reporter lines can lead to data misinterpretation. Here, we identify a consistent, driver-independent tdTomato signal in the adult proventriculus from the widely used lexAop-tdTomato.nls reporter line. This signal was observed across multiple lexA driver combinations and was directly detectable in lexAop-tdTomato.nls responder-alone adult proventriculi lacking any lexA driver and without antibody staining. In contrast, no comparable native red fluorescence was detected in larval proventriculi under the same no-antibody imaging condition. Mouse and rabbit anti-RFP immunostaining further supported the presence of proventriculus-associated tdTomato/RFP antigen in adult responder-alone animals. In larval responder-alone proventriculi, antibody-amplified staining was antibody-source-dependent: a detectable signal was observed only with rabbit anti-RFP, whereas mouse and rat anti-RFP produced no reliable detectable signal under the same staining condition. A driver-matched comparison using lexAop-RFP.nls did not reproduce the proventricular signal, arguing against detectable ectopic activity of the tested lexA driver in this tissue. However, because lexAop-tdTomato.nls and lexAop-RFP.nls differ in reporter/transgene architecture and possibly genomic insertion context, the underlying cause cannot be assigned specifically to the lexAop sequence. Our findings highlight the necessity of including driver-negative and no-antibody controls when using this reporter line in adult Drosophila proventriculus and gut studies.

Gene model for the ortholog of tango (tgo) in the Sep. 2015 (UC Berkeley ASM127793v1/DbusGB1) Genome Assembly (GenBank Accession: GCA_001277935.1) of Drosophila busckii. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

BackgroundThe fission yeast stress-activated protein kinase (SAPK) pathway includes a conserved mitogen-activated protein (MAP) kinase cascade that regulates multiple cellular processes and is activated by several types of external stress. Understanding how Sty1, the MAP kinase in the SAPK pathway, controls these processes is complicated by the fact that different stressors can have stressor-specific effects that may be difficult to separate from the effects of Sty1 activation itself. Moreover, upon stress, Sty1 activation is usually short-lived. Previously, we developed a fission yeast strain, SISA, in which Sty1 kinase activity can be switched on in a sustained manner in the absence of external stress. This required combining multiple mutations in the SAPK pathway, including an analog-sensitive version of Sty1. When SISA cells are grown in the presence of analog-sensitive kinase inhibitors, Sty1 is inhibited, but when inhibitor is removed, Sty1 becomes hyperactive. While this strain was useful, it had several limitations. ResultsHere we describe and validate a more rationally-designed strain, SISA4, that retains the features of the original SISA strain while overcoming its limitations. SISA4 is more stable genetically than SISA, easier to use in genetic crosses, and easy to identify by phenotype or genotyping. We show that analog-sensitive kinase inhibitors 4-Amino-1-tert-butyl-3-(1-naphthylmethyl)pyrazolo[3,4-d]pyrimidine (1-NM-PP1) and 4-Amino-1-tert-butyl-3-(3-bromobenzyl)pyrazolo[3,4-d]pyrimidine (3-BrB-PP1) are equally potent for inhibiting analog-sensitive Sty1 in vivo, and we determine optimal inhibitor concentrations for converting SISA4 cells from a Sty1-inhibited state to a Sty1-hyperactive state. We also find that both 1-NM-PP1 and 3-BrB-PP1 have measurable off-target effects in wild-type cells, although these are modest and generally do not affect interpretation of experiments. Finally, using SISA4, we show that the Sty1-activated transcription factor Atf1 plays an unexpected role in maintaining cell-polarity disruption after Sty1 hyperactivation. ConclusionsSISA4 will be useful for investigating how SAPK pathway activation regulates diverse cellular processes.

Gene model for the ortholog of raptor in the D. grimshawi May 2011 (Agencourt dgri_caf1/DgriCAF1) Genome Assembly (GenBank Accession: GCA_000005155.1) of Drosophila grimshawi. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of raptor in the May 2011 (Agencourt Dere_CAF1/DereCAF1) Genome Assembly (GenBank Accession: GCA_000005135.1) of Drosophila erecta. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

We present Nemo 2.4, an advanced forward-time individual-based simulation framework designed to model the complex eco-evolutionary dynamics and genetic basis of quantitative traits. This tool addresses current challenges in evolutionary quantitative genetics by providing unprecedented flexibility and computational efficiency. Nemo 2.4's modular architecture allows researchers to design custom life cycles by combining specialized Life Cycle Event (LCE) modules, from reproduction and dispersal to selection, crossing, and phenotype expression. The software supports diverse population models, including both Wright-Fisher (WF) and non-WF dynamics, spatially explicit models, and varying demography. Nemo 2.4 handles a wide range of genetic architectures, including both multi-allelic Quantitative Trait Loci (QTL) for general trait studies, and dense di-allelic Quantitative Trait Nucleotides (QTN) implemented with highly optimized bit-wise data structures. Crucially, it allows the simulation of QTNs on comprehensive genetic maps that incorporate other genetic elements, providing genomic-scale resolution. Key biological complexities are integrated natively: the model accommodates modular pleiotropy, dominance, and pairwise epistasis across multiple traits, facilitating the study of complex genotype-phenotype mappings. Furthermore, Nemo 2.4 models phenotypic plasticity through reaction norms and incorporates underlying liability thresholds, enabling the simulation of environmental influences on trait evolution with various forms of selection (e.g., Gaussian, linear, truncation). Due to its compiled design and memory-efficient data representations for large numbers of loci, Nemo provides a robust platform for running high-throughput simulations critical for testing theoretical predictions in polygenic adaptation and understanding evolutionary responses to changing environments.

Developmental processes depend on carefully coordinated gene expression. Expression is modulated by the binding of transcription factors (TFs) to cis-regulatory elements (CREs), like enhancers and promoters. Many computational and experimental approaches have been developed to find CREs, particularly enhancers, in the genome, each with strengths and caveats. Given the increasing availability of ATAC-seq data and methods to find TF binding therein, we hypothesized that we could use TF footprinting tools to find clusters of TF binding events within accessible chromatin that may act as CREs. Using Drosophila anterior-posterior patterning network as a test bed, we used a digital genomic footprinting tool (DGT), TOBIAS, on previously published early embryo ATAC-seq data to characterize the TF footprint landscape of 16 TFs essential for embryonic patterning. Even in this system, with its extensive enhancer annotation, most footprinted TF binding sites lie outside of known enhancers, with intergenic and intronic regions hosting the highest TF footprint count, albeit at low density. To find potential novel enhancers, we identified high-density TF footprint clusters that are highly conserved and overlap with active enhancer histone mark signals. Five high confidence candidates were selected for reporter assay validation and all five were found to drive spatially patterned expression in the embryo. This study shows that even in a highly characterized system, the analysis of footprinted TF binding sites in ATAC-seq data can uncover new regulatory regions and suggests this approach may be helpful in using existing ATAC-seq data to find novel CREs. ARTICLE SUMMARYGiven the increasing availability of ATAC-seq datasets, workflows to exploit the data to uncover new cis-regulatory elements (CREs), including enhancers, are valuable. Using early anterior-posterior patterning in the Drosophila embryo as a test case, we find that previously published transcription factor footprinting tools and ATAC-seq data can be analyzed to yield new candidate CREs. Experimental validation confirms the activity of selected candidate CREs, suggesting that existing data can be analyzed to find novel regulatory elements.

Many biological characteristics arise by interactions between more than one biological organism or unit. Fertilization success in sexually reproducing species represents such an extended phenotype where both mates are required to be fertile for a successful outcome. Consequently, predictive models should account for the joint nature of reproductive performance while offering interpretable estimates for individual mate contributions. Recent advances in genomics and machine learning (ML) provide standardized, high-dimensional genetic information on one hand and computational tools capable of modeling complex biological systems on the other. Here, we construct and evaluate two-tower (TT) machine learning architectures for genomic prediction of binary reproductive outcomes and recovery of sex-specific fertility liabilities. Simulated datasets, generated under a range of genetic architectures, were utilized to compare multilayer perceptron (TT-MLP), convolutional neural network (TT-CNN), and L1-regularized linear (TT-LASSO) two-tower models. Simulation scenarios varied sex-specific heritabilities, genetic correlations, infertility prevalence, mating structure, and sex-specific infertility rates. Models were evaluated with regard to their ability to predict reproductive success at pair level and also recover true underlying genetic values for male and female fertility. Prediction accuracy increased with the underlying heritable component as expected, while sex-specific tower-scores successfully recovered latent fertility liabilities despite models being trained only on observed joint outcomes. TT-LASSO achieved the highest overall classification performance, whereas TT-MLP provided more balanced and consistent recovery of sex-specific genetic values across scenarios. An additional simulation, incorporating genotype-dependent mate compatibility demonstrated advantages of fully-connected neural networks for capturing non-additive interactions. These results indicate that two-tower frameworks provide a powerful approach for modeling reproductive traits, enabling simultaneous prediction of aggregate reproductive outcomes and sex-specific fertility liabilities from genotypic information.

Gene model for the ortholog of Density regulated protein (DENR) in the Apr. 2013 (BCM-HGSC/Deug_2.0) (DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2) of D. eugracilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

The kynurenine pathway metabolizes tryptophan into 3-hydroxykynurenine (3-HK), a precursor for ommochrome eye pigments synthesized via the cardinal (cd) gene in mosquitoes. While cd disruption was presumed neutral, we observed fitness costs in Anopheles stephensi knock-in but not knock-out cd mutants. Here we investigated this anomaly further by assessing survival, fecundity, and midgut integrity across multiple cd mutant lines. Heterozygous knock-in lines, expressing a fluorescent marker and guide RNA for CRISPR/Cas9, exhibited reduced survival post-blood feeding, larva-to-adult survival deficits, and midgut barrier dysfunction, whereas knock-outs showed no such costs. Oral supplementation with xanthurenic acid partially rescued knock-in mortality, implicating oxidative stress linked to 3-HK metabolism. Expression analyses suggest transgene insertion effects, rather than cd disruption, underlie these fitness costs. These findings highlight the importance of evaluating insertional effects in gene drive target selection and support cd as a viable target for genetic control strategies in An. stephensi.

BackgroundStructural variants (SVs) are large insertions or deletions of DNA sequences. While less numerous than single nucleotide polymorphisms, SVs often account for a greater proportion of nucleotide differences between genomes. Their size and frequent association with repetitive sequences has historically hindered their detection, which has limited the ability to associate this variation with molecular and phenotypic trait variation. While some SVs have been linked to observable traits, it remains unclear whether such effects are rare or broadly distributed across the genome. ResultsTo test for genome-wide relationships between SVs and gene expression, we analyzed genome assemblies and transcriptomic data from 10 tissues across 26 diverse maize inbred lines. We identified SVs amongst these lines and examined variants located within the 1kb promoter region upstream of genes. Thousands of genes showed expression differences associated with promoter SVs, often in a tissue-specific manner. One common feature of these SVs was the presence of transposable element sequences. LTR retrotransposons were enriched amongst promoter SVs associated with differential expression and often reduced expression of the nearby gene. Despite widespread expression changes, we found no enrichment for specific biological functions or pathways among affected genes. ConclusionsOur findings indicate that extant TE-mediated promoter SVs play a significant role in shaping gene expression patterns across the maize genome. However, their phenotypic effects appear limited or context-dependent, suggesting that many variants may have minimal impact outside specific developmental stages or environmental conditions.

Genetic-background studies require defined perturbations that can be crossed reproducibly into many recipient backgrounds. We generated a Drosophila dilp2GS-rpr donor line for adult-inducible ablation of insulin-producing cells (IPCs), which secrete insulin-like peptides and provide a tractable model of insulin-deficient metabolic physiology. This line carries dilp2-GeneSwitch-GAL4 and UAS-reaper in cis on the same second chromosome homolog over a balancer. PCR genotyping and sequencing confirmed both transgenic elements in the candidate recombinant line. RU486 induction reduced dilp2 mRNA expression, supporting partial IPC ablation. Treatment-duration testing identified 8 days of RU486 as sufficient to increase whole-body glucose in the dilp2GS-rpr line but not in the background-matched control; food intake did not differ between RU486- and vehicle-treated flies. Across metabolic assays, whole-body glucose showed the clearest RU486- and line-dependent phenotype. This validated dilp2GS-rpr line enables testing how recipient genetic backgrounds modify inducible IPC/DILP metabolic phenotypes and provides a framework for similar linked donor-line resources.

Gene duplication typically fails to confer a selective advantage to an organism, prompting their removal from a population. In the rare instance that duplication either does not incur a fitness cost or it enhances fitness, gene families can form through repeating the duplication process. While the function of gene duplicates has been studied in detail, little work has explored how repeated duplication impacts paralog redundancy and may restrict the emergence of new paralogs or novel function. Here, we constructed a panel of single deletion mutants for each of the 14 members of the <em>Candida albicans</em> telomere-associated (<em>TLO</em>) gene family to test the redundancy in molecular and biological function among paralogs from a lineage-specific expansion. Tlo proteins function as interchangeable subunits of the Mediator transcriptional regulatory complex and have the potential to alter gene expression and an array of cellular responses. Redundancy was the most common outcome, being observed for approximately 80% of the phenotypic assays in strains lacking single <em>TLO </em>genes. However, mutants for all 14 paralogs displayed non-redundant functions in phenotypes ranging from carbon utilization to <em>in vivo</em> virulence. Analysis of gene expression in single <em>TLO </em>mutants found similar trends in redundancy, and loss of single <em>TLO</em>s disproportionately affected genes involved in filamentation, adhesion, redox reactions, and transporter activity at the cell surface. Importantly, sequence divergence between paralogs positively correlated with the frequency of altered phenotypes in single <em>TLO </em>mutants, indicating the acquisition of non-redundant function with increased evolutionary distance. Double mutants lacking two <em>TLO</em> genes produced both positive and negative synergistic phenotypes, suggesting that crosstalk or coordinated regulation is common among paralogs. Together, this study demonstrates that recently emergent paralogs acquire non-redundant functions despite often retaining redundancy with other gene family members to form a highly interconnected functional network.

Animals fight for resources to obtain fitness benefits; most contests are intrasexual, and males tend to fight more than females. Although the genetic basis of male aggression is well studied, we know little about the genetic variation of female aggression. Female aggression varies with reproductive status and is potentially influenced not only by her genotype, but also by the genotype of her mate. Here we measured both male and female aggression in a set of Drosophila melanogaster inbred lines by competing each line against a standard competitor. Aggression varied among lines for both sexes, but male and female aggression were not correlated. Female aggression for many lines increased with mating, as expected, but not all lines changed aggression. However, when females were mated to males of different lines, male genotype did not affect the post-mating change in aggression, suggesting that ejaculate-mediated effects do not vary across these lines. The aggression level of the standard opponent was positively correlated with that of focal individuals indicating that individuals modulate their behavior according to the genotype of their opponent.

Homing endonuclease genes spread by cleaving homologous chromosomes that lack the endonuclease cassette, after which repair from the endonuclease-containing chromosome converts the cut allele into a copy of the drive allele. This mechanism has provided a conceptual foundation for synthetic gene drive systems, including CRISPR-based drives, that represent promising strategies for the genetic control of insect pests. However gene drive performance depends critically on the repair pathways available in the germline of the target organism. Here, we report a set of transgenic assays originally developed as part of an attempt to establish gene targeting in the malaria mosquito Anopheles gambiae using an in vivo-generated linear targeting molecule. Although the intended FLP-mediated excision step was not achieved in the mosquito germline, analysis of the component strains revealed efficient germline activity of the rare-cutting homing endonuclease I-SceI and a striking bias towards homology-based repair of I-SceI-induced double-strand breaks. Across reporter and donor configurations, cleavage outcomes were dominated by single-strand annealing, microhomology-mediated repair, synthesis-dependent strand annealing and gene conversion-like events, with comparatively limited evidence for classical non-homologous end joining. In reciprocal crosses designed to distinguish gene conversion from gamete loss, I-SceI cleavage also produced inheritance distortion consistent with both conversion of the cleaved allele and reduced recovery of gametes carrying extensively damaged donor alleles. These findings indicate that the An. gambiae germline can strongly favour homology-dependent repair following homing endonuclease cleavage and that cleavage can also generate meiotic drive-like distortion through selective loss of damaged gametes. The results have direct relevance for the design and interpretation of homing endonuclease and CRISPR-based gene drives in malaria mosquitoes, where the balance between homology-directed repair, end joining and gamete viability will determine drive efficiency, resistance formation and transmission bias. Author summaryGene drives depend on a simple but demanding principle: a nuclease cuts one chromosome, and the cell repairs the break using the homologous chromosome as a template, copying the drive element in the process. Before CRISPR, this type of system was explored using naturally occurring homing endonucleases such as I-SceI. We attempted to develop a gene targeting system in Anopheles gambiae based on the Rong and Golic strategy, in which FLP recombinase would excise a donor molecule and I-SceI would linearise it to stimulate recombination. The full knockout technology did not work because FLP-mediated excision was not detected in the mosquito germline. However, the component tests revealed something more broadly important: I-SceI-induced breaks were repaired predominantly through homology-based pathways rather than simple end joining. We also observed inheritance distortion consistent with both gene conversion and loss of damaged gametes. These results help explain why homing-based systems can work in mosquitoes, while also highlighting why repair pathway choice and gamete viability need to be measured directly in any new drive configuration.

The P-element, a DNA transposon, has independently invaded two Drosophila species, accompanied by rapid evolution of suppression. In the germline, suppression is mediated primarily by maternally expressed piRNAs, a class of regulatory small RNAs associated with PIWI proteins. The offspring of females that lack P-element-specific piRNAs and males that contain P-elements suffer a syndrome of deleterious phenotypes, including sterility, genome rearrangements, gonadal atrophy, and mutations, while the offspring of the reciprocal cross are normal. These effects, collectively termed hybrid dysgenesis, have been investigated primarily in female D. melanogaster. Here, we study hybrid dysgenesis in male D. simulans. Using an attached-X chromosome stock, we generated genetically identical F1 males that differed only in maternal suppression of the P-element. Using targeted sequencing of P-element breakpoints, we show that P-element transposition is elevated in dysgenic males and confirm a preference for insertion near origins of replication. Using transcriptomics, we show that dysgenic males have elevated P-element expression and reduced splicing suppression, with patterns of gene expression suggesting the loss of mature sperm cells. Fertility assays show higher rates of male sterility but otherwise modest effects on fertility. In conjunction with the transcriptomic data, small RNA sequencing confirms that the piRNA pathway functions in testes. Our results suggest that the P-element may spread more readily through males than females, as transposition rates are similar while fertility defects are less severe in males.

Haplotype-based breeding strategies have emerged as promising approaches to maximize long-term genetic gain by identifying complementary parental combinations while maintaining genetic diversity. However, these methods typically require phased genotypes and more intensive workflow pipelines and skillsets. We developed a novel local genomic estimated breeding value (localGEBV) fitness function with similar intent to the optimal haplotype stacking (OHS) framework fitness function and implemented both in the novel R package, HapSelect. Our aim was to evaluate whether phased haplotypes provide additional benefit over the more easily available dosage-based unphased genotypes in highly inbred crops. A subset of bread wheat nested association mapping (NAM) population comprising 444 lines genotyped with 6,054 DArT-Seq markers was analysed. Marker effects were estimated using rrBLUP, localGEBV and haplotype effects were calculated across linkage disequilibrium-defined haploblocks, and genetic algorithms (GA) were used to identify optimal sets of 30 founders using either a localGEBV derived fitness function with unphased, dosage inputs or the OHS fitness function with phased inputs. Selected parental sets were compared with conventional truncation selection (TS) through 150 generations of forward simulation. The OHS fitness function achieved a marginally greater optimized ultimate GEBV than the localGEBV fitness function during GA optimization, with only 18 of the 30 selected founders overlapped between the two methods. Despite these differences, forward simulations demonstrated nearly identical long-term genetic gain for localGEBV and OHS-selected founders, with both approaches outperforming conventional truncation selection by maintaining greater genetic diversity and delaying the genetic plateau. The minimal difference between localGEBV and OHS is likely attributable to the high homozygosity of the population, where localGEBV and haplotype effects are nearly confounded. These results demonstrate that dosage-based localGEBV provides a practical alternative to phased haplotype approaches for parent selection in inbred crops, substantially simplifying genomic workflows while maintaining long-term breeding performance. Future work should evaluate these methods in more diverse inbred populations and outbred species, where great haplotypic diversity may increase the advantage of true haplotype-based optimizations.

Multi-environment trials are central to cultivar evaluation because they reveal how candidate cultivars perform across locations, years, management conditions, and stress environments. The resulting yield matrix is a rich source of data on genotype-by-environment interaction, and a wide literature on estimation, decomposition, visualisation, and prediction of yield potential and stability has flourished. However the ultimate question of which cultivar to recommend on the basis of such a matrix is often left implicit. The question is far from trivial, and in this paper we formulate cultivar recommendation as an axiomatic ranking problem. This framework is rich enough to encompass the existing literature on stability indices, as well as any other deterministic ranking procedure. We show that many commonly used stability-based procedures can violate minimal criteria of efficiency or consistency. The result of such violations is that a cultivar with uniformly high yield could be ranked below a cultivar with uniformly low yield, or the relative ranks of two cultivars could depend on whether or not a third cultivar is present in the matrix. Our results prove that under a small number of such criteria the space of admissible rules collapses to the family of power means and their limiting cases. If we further wish to allow multiplication normalisation of yield, we are left with the geometric mean as the unique solution.

UNC-6/Netrin is a conserved regulator of dorsal-ventral axon and cell migrations. UNC-6 is composed of a Laminin N-terminal domain (LN), three epidermal growth factor repeats (EGF), and a Netrin C terminal domain (NC). Here, we identified missense mutations in distinct UNC-6 domains and assessed their roles in dorsal VD/DD motor axon guidance and ventral AVM axon guidance. A missense mutation in a conserved residue of the LN domain (G289D) resulted in dorsal and ventral axon guidance defects similar to unc-6 null. A distinct missense mutation in the LN domain (S120F) was hypomorphic and strongly perturbed ventral AVM axon guidance with minimal effects on dorsal VD/DD axon guidance, showing that S120F is predominantly required for ventral guidance. Missense mutations altering conserved cysteine residues involved in di-sulfide bonding in the EGF domains were analyzed. EGF1(C321G) caused both ventral and dorsal axon guidance defects albeit weaker than unc-6 null, indicating that EGF1 is required for both. EGF2(C347Y) strongly affected dorsal VD/DD axon guidance similar to unc-6 null, with weaker perturbation of ventral AVM axon guidance. Previous results revealed that EGF3(C410Y) specifically disrupted dorsal axon guidance, a result that we confirmed. Our studies using missense mutations in the endogenous unc-6 locus complement previous structure-function studies using transgenic expression, and identify domains specifically required for ventral AVM guidance (S120Y in the LN domain) and dorsal VD/DD axon guidance (C410Y in EGF3). The crystal structure of UNC-6 indicates conserved N-linked glycosylation at N114 and N128. Mutation of these sites in UNC-6 had no effect on dorsal ventral axon guidance, showing that they do not play a major role. However, the N114 and N128 mutations interacted genetically with unc-40 and unc-5 mutations, indicating that these glycosylation sites indeed have a role in UNC-6 signaling. Our results will inform studies on how these distinct UNC-6 domains interact with guidance receptors (e.g. UNC-40/DCC and UNC-5) and other extracellular molecules to mediate dorsal-ventral axon guidance.

If an individuals niche is determined by its genotype, then competition for limiting resources should be most intense among individuals of the same genotype. Theory predicts this will act to maintain genetic variation, but whether this mechanism operates under natural conditions is unclear. Using long-term observations of a population of free-living Trinidadian guppies, we asked (i) whether competition was strongest between kin, and (ii) whether this process maintained genetic variation. Competition between kin was 1.5-1.8 times stronger than that between non-kin. This contributed to balancing selection: after [~]10 generations, variation was 29% higher than expected under drift. Our results show that relatedness can play a major role in structuring ecological competition, with broader consequences for theories of inclusive fitness. One-sentence summaryHeritable variation is maintained due to resource competition being more intense among kin.